Statistics of the Mesoscopic Field 



in 
o 
o 

(N 



A.A. Chabanov^ 2 -0 and A.Z. Genack 1 

1 Department of Physics, Queens College of the City University of New York, Flushing, NY 11367, USA 
d CEMS, University of Minnesota, Minneapolis, MN 55455, 
(Dated: 14 February, 2005) 



USA 



We find in measurements of microwave transmission through quasi-lD dielectric samples for both 
diffusive and localized waves that the field normalized by the square root of the spatially averaged 
flux in a given sample configuration is a Gaussian random process with position, polarization, 
frequency, and time. As a result, the probability distribution of the field in the random ensemble 
is a mixture of Gaussian functions weighted by the distribution of total transmission, while its 
correlation function is a product of correlators of the Gaussian field and the square root of the total 
transmission. 
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Enhanced fluctuations and correlation of the local and 
spatially averaged electronic and electromagnetic flux in 
mesoscopic samples have been extensively studied in the 
last two decades Lil L3- Lit ■ Yet the impact of weak and 
strong localization upon the statistics of the field itself 
has not been investigated. The field at a point in a 
multiple-scattering medium is the superposition of par- 
tial waves associated with all wave trajectories or Feyn- 
man paths reaching that point. The changing distribu- 
tion of paths with position produces a speckled intensity 
pattern which fluctuates on a scale of the wavelength, A. 
When a large number of statistically independent terms 
contribute to the field and the sample is far from the 
localization threshold so that long-range correlation is 
not substantial, the distribution of in- and out-of-phase 
components of the field is Gaussian JiJ. This leads to 
a negative exponential distribution for the intensity of 
a polarized component of the field, which is widely ob- 
served in diffusive media. 

In mesoscopic samples, however, the crossing of Feyn- 
man paths gives rise to correlation in the intensity at 
distant points. This leads to enhanced fluctuations of 
intensity 0, El HI Q l , total transmission 0, HH Ull ] , and 
conductance j 1 2L Il3| , with increasing variances as the lo- 
calization threshold is approached. Since the statistics of 
transport in mesoscopic samples reflects the statistics of 
the underlying field, it is of interest to investigate field 
statistics in samples with increasing intensity correlation. 

In this Letter, we investigate the structure of the prob- 
ability distribution and correlation of the mesoscopic field 
in both the frequency and time domain for diffusive and 
localized waves. We find that the distribution of the field 
in a given realization of a random quasi- ID sample and 
in the subset of realizations with the same value of to- 
tal transmission is Gaussian, even for strongly localized 
waves. The field distribution in the random ensemble 
is therefore a mixture of Gaussian functions weighted by 
the distribution of total transmission. Thus the statistics 
of the total transmission can be obtained from the statis- 
tics of the field, as well as the reverse. When the field is 



normalized by the square root of the spatially averaged 
flux in each configuration, it becomes a Gaussian random 
process with position, polarization, frequency, and time. 
The statistical independence of the normalized Gaussian 
field and square root of the total transmission allows the 
field correlation function to be written as a product of 
their respective correlation functions. The field corre- 
lation function with frequency or time shift exhibits the 
impact of localization through the corresponding correla- 
tion function of the square root of the total transmission, 
while the field correlation function with displacement and 
polarization rotation has a simple universal form in both 
the frequency and time domain for diffusive and localized 
waves. 

To examine field statistics, we measure microwave 
spectra of the field transmission coefficient in ensembles 
of random quasi-lD samples of alumina spheres with the 
use of a vector network analyzer. Alumina spheres with 
diameter 0.95 cm and refractive index 3.14 are embedded 
in Styrofoam spheres of refractive index 1.04 to produce 
an alumina volume fraction of 0.068. The spheres are 
contained in a copper tube of diameter 7.3 cm and length 
61 cm with reflecting sidewalls and open ends. Spectra 
are taken in 10 4 sample configurations produced by rotat- 
ing the sample tube, in steps of 0.3 MHz in the frequency 
intervals 14.7-15.7 GHz (measurement A) and 9.95-10.15 
GHz (measurement B) , in which waves are diffusive and 
localized, respectively |l4j. The frequency intervals are 
sufficiently narrow that propagation parameters change 
little within each interval. 

The probability distributions of the transmitted inten- 
sity and total transmission normalized by their ensemble 
averages, s ab =\t ab \ 2 / (\t ab \ 2 ) and s a =S h |i ah | 2 /(E fc |i afc | 2 ), 
respectively, where t ab is the transmission coefficient for 
incident field a and outgoing field b, are determined 
by the variance of the normalized total transmission, 
var(s a ) 0, [n|. Var(s a b) can be found from the rela- 
tion, var(s a £,)=l+2var(s a ) Q. In the limit of Gaussian 
field statistics for the random ensemble, var(s a b)=l 0]. 
In measurements A and B, var(s a &)=1.18 and 6.18, re- 
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spectively. In order to study field statistics in even more 
strongly correlated samples, we examine the statistics of 
pulsed transmission of localized waves for a long delay 
from an exciting pulse, since the degree of correlation in- 
creases with time delay 0, . The response to a pulse 
with a Gaussian temporal envelope of width a t =160 ns is 
obtained from the Fourier transform of the product of the 
field transmission spectra of measurement B and a Gaus- 
sian spectral function of bandwidth er=lMHz w 0.6(5^, 
where 8v is the field correlation frequency ^3] ■ The field 
statistics are examined for waves delayed by 740 ns from 
the center of exciting pulse. In this case (measurement 
C), var[s o6 (f)]=20.1. 

We first consider the probability distribution of the 
field transmission coefficient normalized to the square 
root of the ensemble-averaged intensity, E=t a t, j\J (|i o| 2 )- 
We find in all cases that the statistics of the real and 
imaginary parts of the field, r=Kc[E] and i=Im[E], re- 
spectively, are the same and that r and i have zero 
mean and vanishing cross correlation. The distributions 
P(a), where a—r,i, are shown by the solid curves in 
Fig. 1. Increasing deviations from a Gaussian distribu- 
tion, P(a)=-^ exp(— a 2 ), obtained for the model of the 
field as a random phasor sum (dashed curve) , are seen 
for the distributions in measurements with larger values 
of var(s ab ). 

Though departures from a Gaussian distribution in en- 
sembles of mesoscopic samples arise from extended cor- 
relation generated by the crossing of wave trajectories, 
we expect that the field distribution in a given quasi- ID 
sample configuration is Gaussian. Once wave trajecto- 
ries cross in the medium, intensity fluctuations propa- 
gate throughout the sample, leading to correlation in the 
field amplitude on the output surface at points separated 
by considerably more than the field correlation length of 
A/2 This gives rise to enhanced fluctuations of 
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FIG. 1: Field distributions P(a) of measurements A, B, and 
C (solid). The dashed curve is the Gaussian distribution, 
P(a) = -j- exp(— a 2 ), of the model of the field as a random 
phasor sum 0. 



the total transmission from configuration to configura- 
tion. However, in a given configuration the field at each 
point on the output surface is the sum of a large number 
of partial waves with identically distributed amplitudes 
and phases, which is due to perfect mixing of modes in 
a quasi-lD sample. When the number of field coherence 
areas on the output surface is large, the field in a given 
configuration and in the subset of configurations with a 
specific value of s a is a Gaussian random variable. We 
therefore expect that the joint distribution of the real 
and imaginary parts of the field in the random ensemble, 
P(r, i), will be a mixture of Gaussian distributions [2(J of 
zero means and variances s a /2, with mixing proportions 
given by the total transmission distribution P(s a ), 



f°° 1 

P(r,i) = / ds a P(s a ) exp 

Jo n s a 



r 2 + i 2 



(1) 



Since s ab = \E\ 2 = r 2 + i 2 , the distribution of normal- 
ized intensity, P(s a b), is a mixture of negative exponen- 
tial distributions, 



P(sab) = [ ds a P(s a )— exp [-— J 

JO s a \ S a / 



(2) 



This result is consistent with random-matrix calculations 
carried out by Kogan and Kaveh 0. 

The field of Eq. (1) may be expressed as the product 
of two statistically independent random variables, E= 
E'x^/s^, where E'=Ej^fs2 is a Gaussian field with dis- 
tribution P{r' exp[-(r' 2 +i' 2 )], where r'=Re[E'] 
and i'—lm[E']. This can be demonstrated by obtaining 
the distribution P(r,i) of Eq. (I) from the joint distri- 
bution P(r',i';s a ). Since P(r', i'; s a )=P{r' , i')xP(s a ), 
the probability $ e {X, Y) that the real and imaginary 
parts of E assume values less than or equal to X 
and Y, respectively, can be expressed as &e(X, Y) = 
J °° $ B /(-^=, -j=)P(s a )ds a - Its derivative is the joint 

probability density function, Pe(X,Y 

Jo 7^ r E>{^ 
with Eq. (1). 

The normalized intensity s ab can also be written as 
a product of two statistically independent variables, 



d 2 <S> E (X : Y) _ 
dXdY 

-¥=)P(s a )ds a , which is in agreement 



where s' 



E 



l\2 



The corre- 



sponding joint distribution is P(s' ab ; s a )= P(s' ab )xP(s a ), 
where P(s^ b )=exp(— s' ab ). Since the n-th moment of the 
product of two statistically independent random vari- 
ables is the product of their n-th moments, the moment 
relation between the field and the total transmission is 



(2fc-l)!! / fc\ 
2*= \ s a/ 



0, 



n 
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2k 
2k 
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(3) 



where a'=r',i' and k is an integer, and the moment re- 
lation between the intensity and the total transmission 
is 
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Equation (4) is in agreement with the calculations of 
Ref. and measurements of Ref. [f| . 

To demonstrate that P(r, i) is the mixture of Gaussian 
distributions of Eq. (1), we solve for P(s a ) by utilizing 
the field average (cos(2ay / i)), with z>0. Expanding the 
cosine as a power series in a and then using the moment 
relation of Eq. (3), we obtain the Laplace transform of 
P(s a ),F(z), (cos(2av^)) = (exp(-s Q z))=F(z). The solid 
curves in Fig. 2a are plots of F(z)—(cos(2a^/z)) versus z 
for measurements A, B, and C. These curves are inverted 
to obtain P(s a ), b y u s ing approximate inversion of the 
Laplace transform j2lj,|22|. The distributions P(s a ) are 
shown by the solid curves in Fig. 2b. 

F(z) is also equal to the intensity average 
(Jo(2^/zs a fc)), where Jo (a;) is the Bessel function of 
zero order. This can be shown by expanding the Bessel 
function as a power series in s a b and then using the 
moment relation of Eq. (4). F(z)—(J (2y/zs a b)) for each 
measurement is plotted as the dotted curve in Fig. 2a 
and coincides with F(z)=(cos(2a^/z)). 
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FIG. 2: (a) Laplace transform F(z) of the total transmission 
distribution P(s a ), found as the field and intensity average, 
F(z)={cas(2a^/z)) (solid) and F(z)=(Jo(2 s /zs^)) (dotted), 
respectively, in measurements A, B, and C. The dashed curves 
are F(z) of Eq. (5) [aQ with g replaced by 2/3var(s a ). (b) 
P(s a ) found by inverting the corresponding F(z) in (a). For 
measurement A, F(z) is fit by F(z) = exp[Q(z)], where Q(z) 
is a polynomial of power 3, and the resulting function is in- 
verted using the Weeks method [2lll ; for measurements B and 
C, F(z) is inverted using the Stehfest method (2^. 



F(z) was previously found in diagrammatic 6] and 
random-matrix-theory |jj calculations in the diffusive 
limit <?3>1, in the absence of absorption, to be a func- 
tion only of the dimcnsionlcss conductance g, 

F{z) = exp[-gln 2 ( v /l + z/ 3 + yffig)]. (5) 

A nonperturbative result for F(z) was obtained in Ref. || 
for the case of broken time-reversal symmetry. Theoreti- 
cal expressions for P(s a ) and P(s a b) derived using Eq. (5) 
were found to closely match the measured distributions 
[ol In! H5I ] , even at the localization threshold, g ~ 1, and 
even in the presence of absorption, once g was replaced 
by 2/3var(s a ). Var(s a ) can be found from the relation, 
var(s a )=i[var(s a f,) — l], following fromEq. (4). The plots 
of F(z) of Eq. (5) with g replaced by 2/3var(s a ) with the 
values of var(s a ) of 0.09, 2.59 and 9.55 for measurements 
A, B, and C, respectively, are displayed as the dashed 
curves in Fig. 2a. The theoretical curves are seen to 
deviate from F[z) found as the field and intensity av- 
erage at large z. Deviations also appear in P(s a ) for 
s a < 1 (Fig. 2b); increasing deviations from the approx- 
imate theory of Ref. 0,0 for P(s a ) are observed for the 
distributions with larger values of var(s a ). 

We now consider field and intensity correlation versus 
shifts in the position and polarization of the detected 
wave, p, and in the incident frequency or positions of 
scatterers within the medium, q. By virtue of the sta- 
tionarity and statistical independence of the normalized 
field E'(p,q) and the total transmission s a (q), the field 
correlation function with shifts in p and q, Ap and Aq, 
F E (Ap, Aq)={E(p, q) E*(p+Ap, q+Aq)), can be written 
as the product of two correlation functions, 

F E {A Pl Aq) = F E ,(Ap, Aq) x T^(Aq) , (6) 

where F E > ( Ap, Aq)=(E'(p, q) E'*(p+Ap, q+Aq)) and 
T^(Aq)=(y/s a (q) s a (q + Aq)). Similarly, the inten- 
sity correlation function, T Sab (Ap, Aq)—(s a b(p, q) s a b{p+ 
Ap, q+Aq)), can be written as 

T Sab (Ap, Aq) = T <b (Ap, Aq) x T Sa (Aq) , (7) 

where IV (Ap, Aq)=(s' ab (p, q) s' ab (p+Ap, q+Aq)) and 
r Sa (Aq)=(s a (q) s a (q+Aq)). The Siegert relation for the 
Gaussian random process E'(p,q) gives T s i (Ap, Aq) — 
l+F'(Ap,Aq), where F'(Ap, Aq)=\F E , (Ap"Aq)\ 2 Q. 

For a shift in position or polarization, Ap, with 
Aq=0, r v ^ r (Ag)=(s a )=l and T Sa (Aq) = (s 2 a ) . This gives 
F E (Ap)=F E ,(Ap) and T Sab (Ap) = ( S 2 a )[l+F'(Ap)]. For a 
shift in position, Ap= Ar, the functional form of F E i ( Ar) 
is that predicted by coherence theor y 115 , as found 
in microwave fljl |23| and numerical |24| studies of field 
correlation. The cumulant intensity correlation func- 
tion with displacement is given by C(Ar)=r Sa6 (Ar)— 
l=F 1 '(Ar)+var(s Q )[l+F'(Ar)]. Identifying the degree 
of nonlocal correlation, the value of C when 
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F'(Ar)=0, we find k is precisely equal to var(s a ) and 
thus C(Ar)=F'(Ar)+K[l+F'(Ar)] [lHH. Similar re- 
sults were found for correlation with rotation of the polar- 
ization of the detected field, Ap=A0 ^5|. For any value 
of var(s a ), Fe'(A9)=cos(A9). The cumulant intensity 
correlation function is C(A6)=F'(A9)+k[1+F'(A9)}. 

When the total transmission s a varies with the inci- 
dent frequency v, or with time r as the internal struc- 
ture of the sample changes, mesoscopic correlation is no 
longer represented by a single parameter. For exam- 
ple, the field and intensity correlation functions with fre- 
quency shift become Fe{Av)=Fe'{Av)xT^-^(Av) and 
T Sab {Av)=[l+F'{Av)}xY Sa {Av), respectively The field 
correlation function Fe(Av) is of particular interest be- 
cause it is the Fourier transform of the time-of-flight dis- 
tribution, P(t), where t is the time delay following a short 
pulse [17|. Since Fe{Av) is the product of two functions, 
P(t) is the convolution of their two Fourier transforms. In 
addition, the correlation function r /^(Af) can be writ- 

ten as T^Au^iy/^f+ST^Au), with ST^Av) 
vanishing at large Av. We therefore expect that P(t) is 
the sum of two terms. The first term is associated with 
spectral correlation of the Gaussian field E', and the sec- 
ond is due to correlation of the square root of the total 
transmission. 

In geometries other than quasi- ID, such as a slab, in 
which transverse modes are not completely mixed, we 
may replace the total transmission s a with a variable 
s Q (r), which is the average intensity within a region on 
the output surface over which field statistics is nearly 
Gaussian. The field and intensity distributions are then 
still given by Eq. (1) and (2), respectively, and their spa- 
tial correlation functions are the products of two terms. 
The first term is a narrow function of Ar, describing cor- 
relation within the speckle pattern, while the second term 
is a broad function, reflecting correlation of the more 
slowly varying function s a (r) over the output surface of 
the sample. 

In conclusion, we have found the manner in which 
field statistics arise from fluctuations and correlation of 
the spatially averaged flux in mesoscopic samples. We 
find that the probability distribution of the field in a 
given quasi- ID sample configuration and in the subset of 
configurations of flux s a is Gaussian and that the field 
distribution in the random ensemble is the mixture of 
Gaussian functions weighted by the distribution P(s a ) 
[Eq. (1)]. Thus, non-Gaussian mesoscopic statistics of 
the field arise from the breakdown of ergodicity in which 
the Gaussian field statistics for different points within 
a given configuration differ from the mesoscopic statis- 
tics at a fixed point in a random ensemble. When the 
field is normalized by in each configuration, it be- 
comes a Gaussian random process with position, polar- 
ization, frequency, and time. The field correlation with 
displacement and polarization rotation is independent of 



closeness to the localization threshold or of the degree of 
correlation, n. In contrast, the field correlation function 
with frequency or time shift is written as a product of 
correlators of the Gaussian field and yfsZ- Since the time- 
of-flight distribution for particles is the Fourier transform 
of the field correlation function with frequency shift, the 
increasing suppression of transport with time delay as 
a result of the increasing impact of weak localization is 
associated with mesoscopic fluctuations. 
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